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Abstract 

In the absence of interaction it is well known that the Kogut-Susskind 
regularizations of fermions in the spin and flavor basis are equivalent to each 
other. In this paper we clarify the difference between the two formulations in 
the presence of interaction with gauge fields. We then derive an explicit ex- 
pression of the transfer matrix in the spin basis by a unitary transformation 
on that one in the flavor basis which is known. The essential key ingredient 
is the explicit construction of the fermion Fock space for variables which live 
on blocks. Therefore the transfer matrix generates time translations of two 
lattice units. 



1 Introduction 



The naive discretization of the Dirac equation on the lattice [1] leads to a replica- 
tion of the fermionic states, known as lattice fermion doubling [2]. The doublers 
appear as spurious poles in the fermion propagator at the nonzero corners of the 
Brillouin zone. The Wilson way of removing the doublers is to give them a mass 
which becomes infinite in the continuum limit, at the cost of an explicit breaking 
of chiral invariance on the lattice [2]. 

In the Kogut-Susskind [3-6] lattice formulation for relativistic fermions [7, 
Chap. 4] the doublers are instead interpreted as physical fields by the introduction 
of additional quantum numbers. This has been done in two ways. In the former 
approach, first the fermion field is reduced to a single component per site by a 
procedure called spin diagonalization, and, for this reason, this method is referred 
to as the one in the spin basis. Afterwards spin and flavour degrees of freedom are 
associated to different corners of an elementary hypercube on the lattice [8-10], 
and therefore sometimes fermions in this formulation are said to be staggered. In 
the latter approach [11-13], said in the flavour basis, the additional quantum num- 
bers, called taste, are associated, together with the spin, with blocks corresponding 
to the hypercubes of the spin basis of size twice the lattice spacing. 

In the absence of coupling with gauge fields these forms are changed into one 
another by a linear transformation on the fermion fields, but in the presence of 
gauge fields they are not equivalent, as we shall make clear in the following. Their 
difference is of consequence in the construction of the corresponding transfer ma- 
trices. 

For Kogut-Susskind fermions in the flavour basis a simple operator realization 
of the transfer matrix is known [14]. It has been built in close analogy with the 
case of Wilson fermions [15-19] (see also [20]), the only difference being that it 
performs time translations by one block instead of one lattice spacing. 

The situation is more complex for Kogut-Susskind fermions in the spin ba- 
sis [11,12,21], because all attempts at constructing a positive definite transfer 
matrix that performs time translations by a single lattice spacing failed. The dif- 
ficulty was circumvented by looking at time translations by two lattice spacings. 
Here we meet with a subtlety. We must distinguish whether the Fock space is built 
on one or two time slices. In the first case, the square root of the matrix which 
translates by two lattice spacings is the one that translates by one lattice spacing. 
In the second case, instead, translations by one lattice spacing are not defined at 
all. This seems to be the case with Kogut-Susskind fermions, but the necessary 
construction of the Fock space on blocks, in the spin basis, has not been made 
explicit. 

We became interested in a formulation of the transfer matrix in the spin basis 
in the framework of relativistic field theories of fermions whose partition function is 
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dominated by bosonic composites ['^■^j. This subject became for us more relevant in 
the development of an approach to QCD hadronization (meant as the replacement 
of the QCD degrees of freedom by hadronic ones) that makes use of the operator 
form of the transfer matrix [23-26]. Using Kogut-Susskind fermions, because of 
the lack of a convenient formulation of the transfer matrix, we were able to express 
our results only in the flavour basis. Numerical simulations are, instead, usually 
performed in the spin basis, because they are much faster. We were thus motivated 
to find an operator form of the transfer matrix in this latter basis as well. Since 
apparently in any case we should resign to time translations by one block, we 
decided to get an expression of the transfer matrix in the spin basis by a linear 
transformation from the flavour basis. 

We deem that the question might be of more general interest, and therefore we 
report our results in the present paper. In Sect. 2 we remind for the convenience of 
the reader and in order to establish the notation what is relevant for the following 
about the Kogut-Susskind regularization. We adopt the notations of Montvay and 
Miinster [7] with some minor changes that will be specified. In Sect. 3 we perform 
the transformation of the action from the flavour to the spin basis. Most of the 
results, with some qualification, are well known, but we think this Section is a 
necessary preparation for Sect. 4, in which we perform the transformation of the 
transfer matrix. 



2 Kogut-Susskind fermions 

Let Xf^ be the coordinates of hypercubic lattice sites, 0<x^<L^ — l,0</i<3 
(Montvay and Miinster in [7] use indices from 1 to 4), and the coordinates of 
hypercubic blocks. They are related by 

Xf, = + T]^ (1) 

with < < L'^ — 1, = 2L^, and rj^ = 0, 1 the position vectors within the 
block. The sum over lattice points can be split into the sum over the blocks and 
the sum over the sites within a block, that is 

X y r] 

We denote by tp^ the fermionic fields on the lattice sites, and by g^" the fields on 
the blocks. The latter have Dirac spinor indices 1 < a < 4, in greek letters, and 
taste indices 1 < a < 4, in latin letters. 

It is important to remark that the gauge transformations in the first case act 
at the sites of the basic lattice, in the second at the coordinates of the blocks 

^x^gxi^x, qT^dyqT- (3) 
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While gy is the same transformation for all a; in a given block with coordinate y, 
will in general change also within the same block. 



2.1 The flavour basis 

The gauge link variables on the blocks are denoted by Un{y). Under gauge trans- 
formations they change according to the rule 

U,(y) ^ gy U,{y) . (4) 

The action of the fermion fields in the flavour basis can be written as 

= 2^5]£,([/) (5) 

y 

where the factor 16 keeps into account the volume of the elementary cell when 
using variables defined on the blocks, and the Lagrangian in the fiavour basis is 



the flavour matrices are deflned for = 0, . . . , 3 and /i — 5 hy 

and the other operators are deflned in terms of translations on the blocks 

y' 

and the identity on the blocks 

My ^= 2^ E^V,. /(?/) = /(?/) 



(6) 



(7) 



(8) 



y' 



(9) 



according to 



Vi+) := I {U,T^;^ - 1) , V(-) := i (1 - T(-)t/t) 
^(vW-V!.-0 = ^(C/.rW+T(-)C/t-21L). 



(10) 
(11) 
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We can recognize that the projections of the fermionic field 

q+ = P+q, ql=P-q (12) 

where 

P± = -{l®lTlol^®hh) (13) 

propagate forward/backward in time, and therefore describe particles/antiparticles 
respectively. Accordingly we introduce creation and annihilation operators q\-,q±- 
They are defined at one and the same time, so that in addition to spin and fiavour, 
they depend on the spatial position only, denoted by boldface letters. They satisfy 
canonical anticommutation relations 

{{qi)Z. {q±t\] = \ ^y.y.i^f •"'^ , {(d)yr> (^t)S} = . (14) 

As the factor 1 accounts for the spatial volume of the blocks, the above anticom- 
mutation relations become canonical in the basis in which P± are diagonal. 

The transfer matrix corresponding to the fiavour-Lagrangian (6) in the gauge 
f/o = 1 is [14, 27] 

Tt,t+i = exp (g_ Nt q+ exp(2/i n^) exp (g_ Nt+i q+ ) . (15) 

In the above equation Nt is a matrix which depends on the time of the blocks only 
because it depends on the gauge link variables 

Nt := N[U{t)] , (16) 

and 11 is the chemical potential 

^B = '2'Y,{q\q+-qU-)^ (17) 
y 

that we omitted for simplicity in the Lagrangian. By keeping into account the 
spatial volume factors 

g_ Ntq^ = 64 J] {q.)ANt)yy{q+)y (18) 
y'.y 

3 

Nyy = -2 |m (7o ® 1) ly'y + 5]](7o7fe ® 1) 

k=l 

where 

) = -(1® I±7fc75®t54) (20) 
5 



V' V 



(19) 



and 



Vy^=^^y'y, (^^)?y^=^W (21) 



enter in the definitions of v[.'^'*. 
Notice tliat 

qiNq± = 0. (22) 



2.2 The spin basis 

For tlie sake of later comparison we report the regularization of a Lagrangian in 
the spin basis. The gauge fields on the hypercubic lattice are denoted by m^(x) 
and transform according to 

Mx) gxu^{x)gl^f^. (23) 
The Lagrangian in the spin basis is 

C^{u) := m^^tp^ + 2 Xl"^/^ [i^xUf,{x)^^+fi - i^^+^ul{x)^^] (24) 

At=0 

where the signs axfj, are defined for /i = 0, . . . , 3 by 

a,, := (_i)-o+-+x.-i . (25) 

There is no direct way of identifying forward and backward movers. This is the 
difficulty encountered in the construction of a transfer matrix in operator form for 
this Lagrangian. Indeed, as far as we know, such a construction has been achieved 
only after a reduction of the Lagrangian itself, in which the fermion fields and 
their conjugates live on odd and, respectively, even sites [ ,]. 

At the classical level, however, the fields in the spin and fiavour basis are related 
according to 

C = ^ Er.;-^2,+, (26) 
C = ^E^2,+,rt^.« (27) 

where 

r, := iTiTiTiT ■ (28) 
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The matrices T satisfy the relations 

1 



^tr(rtr,,)=5w (29) 

7 X] ^U/9^^-°« ^ ^baSpa (30) 
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that allow us to invert Eqs. (27) 

7A2,+,=2tr(rtg,) (31) 
'ip2y+ri =2tr (g^r^) . (32) 

We will use these relationships in order to derive an action and a transfer matrix 
in the spin basis from those in the flavour basis. 

3 Transformation of the Lagrangian 

In this Section we express the Lagrangian (6) in the spin basis using the transfor- 
mations (27) 

5^£;(t/):=2^^£,(t/). (33) 

X y 

While in the absence of gauge interaction C'^ coincides with reported in (24), 
we shall see that this does not occur, in general, in the presence of gauge fields. 
The mass term of the action is proportional to 

J^^y'iy = ^2j/+'y (rj,r^) tP2y+v' = ^^^^ ■ ^^^^ 

y y V iT X 

In order to derive the kinetic term we shall use the relations 

^ ^ ^rj-.aa ^Ori^^-^r/^^ r]+ll:a' a ~l~ ^Ir/ii^-^ri^^ri—fi-.a'a (35) 

a 

—fj,:a'a' 

(36) 

a,a 

From the definition (28) soon follow the relation (35) and 

r,7^ = (-l)vo+vi+V2+vs^_iy,^^Y^ (37) 

so that 

r,75 = (-i)*+''i+''^+''^75r, (38) 
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and therefore 

a, a 

(39) 

which together with (35) imphes the relation (36). 
The kinetic term is proportional to 

y M 

-gy(75 ® ^5) [U^,{y) qy+ii + Ul{y- fi) q-y^j, - 2qy\ ] (40) 

that is 

y tJ, rfyT]' a,a' ,a,a' 



which is because of (35) and (36) 

gEEE E '^■iy+V''^l':a'a'(^r,i^[U^,{y)So^ri^,'^^+f,-.a'a'i^2y+2fL+ri 
y II Ti,Tj' a,a',a,a' 

~Uj^{y ~ (i)hr,fX r)-fi:a' a'i^2y-2fi+r, + {—^Q-qJ^-q+fi-.a'a' + hr^J^ r^- jl:a' a')'^2y+r,\ (42) 

and performing the trace on spinor and flavour indices (30) 
2 E E E ^sjz+V^w \U^l{y)5Q'n^5'n','n+^i'ip2y+2fl+r) 

y tM ri,rj' 

-Uliy - P')5lri^5r,',r,-fi1p2y~2fi+ij + (-^OTy^. + ^Iri^^n' ,n-fi)i'2y+ri] (43) 

and performing the sum over 7]' 

9 E E E [ '^0r,^^2y+„+A U^{y) 1p2 

-5lr,,>2|y+,,-A Uliy - A) ^2(y-f.)+ri " ^o^^^ V'2y+r,+A '^2y+n\ ■ (44) 
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y V 



Remark that if we increase the component of a site x we jump on block 
different from that of x if is odd. This is the case when x = 2y + t] + fi and 
rjfj, = 0, but not when x = 2y + rj — jl and = 1. Similarly, if we decrease x^ we 
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jump on a different block only when is even. This is the case when x = 2y+ri—fi 
and 1]^ = 1, but not when x = 2y + ri + fi and = 0. And that, ii x = 2y + r] then 



a. 



= a 



(45) 



Then the kinetic term has the form as that of C^{u') where 




Uf^{y) for X = 2y + T] and rj^j^ = 1 
1 elsewhere 



(46) 



that is the gauge field couples only sites which belong to different blocks. 
In conclusion 



We have the constraint, however, that the fermion fields within a block should all 
transform in the same way under gauge transformations. One might think that 
we could relax this constraint by a different transformation from the spin to the 
flavour basis 



Such a generalization, however, is only apparent, because the curvature for the 
plaquettes with all the vertices within one and the same block vanishes. Indeed, 
such a generalization, as the particular ones chosen for example in [0, Eq. (35)], [27, 
Eq. (56)] is a pure-gauge transformation of (27). 

We conclude that, in the presence of a generic gauge-field configuration, the 
Lagrangian in the spin basis C'^{u) and that in the flavor basis Cq{U) are not 
equivalent. 

The transformed Lagrangian C'^iv!) could also be regarded, in the spirit of the 
previous quoted attempt [J J ], as a modification of C^{u), defined in (24), for which 
a transfer matrix can be constructed. 

The above construction refers to the case of vanishing chemical potential. Its 
inclusion is, however, straightforward [ ]. We only note that, at variance with 
respect to the coupling with gauge fields, the chemical potential can be attached 
to all links in the transformed Lagrangian C'^{u'), provided its value be half the 
one in the flavour basis. 



C'^{u') = mij^ij^ + 2 XI "^/^ - ^^+/,M2(a;)^x] • (47) 



/i=0 




(48) 



V 
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4 Transformation of transfer matrix and coher- 
ent states 

As a first step we must transform creation-anniliilation operators from the flavour 
to tlie spin basis. To tliis end we must determine tlie expressions of tlie fields q± 
in tlie spin basis 

{q+)y = ^r^^2y+r, {ql)y = P-^^Trtij2y+v ■ (49) 

V V 

Using the relation (39), we find 

P+Tfj = 6orio^ri , P-^-q = ^l-qo^ rj (50) 

and similar relations hold for F^. 
We therefore have 

i<l+)y = ^0r,a'^v^2y+^ , ^'i-)y = \^ hr^o^ri'^p2y+r^ ■ (51) 

Next we define the operators corresponding to the ^/'-fields according to 

(^+)?/ = ^ 5Z ^^m'^v'^^y+v ' ^^-^y = \^ hrio'^v'^-iy+v (52) 
and assume that 

{V^5y'+^', V^2y + J = 2 5y/y . (53) 

This is obviously consistent with the second set of equations in (14). Consistency 
with the first set requires that 

1 X ^ ^ 1 

= lSyyS^rPr^' (55) 

where a = 0, 1 respectively when the index of the projector is + or — . The second 
equality follows from the equations 



a' ,a' Tj r) 

that can be proven using (30) and (50). 
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Some comments about our results are in order. We see that the temporal com- 
ponent 7]q of the fields in the spinor basis corresponds to the ± projection of the field 
in the fiavor basis. The 8 Dirac-taste degrees of freedom of particles/antiparticles 
are spread on the 8 sites of the even/odd time slice in the corresponding block. 
In this connection, looking at Eq.(53), tjq can be regarded as a quantum number. 
But this quantum number changes when time increases by one unit in the original 
lattice, so that, unlike the q± projections, the fields ip2y+r] with r]Q respectively 1 
or cannot be identified as forward/backward movers. Changing time we change 
a particle into the hole of an antiparticle. 



4.1 Transfer matrix 

We first transform the baryon number 



(57) 
(58) 



where we re-label the operators with the spatial coordinates 

x = 2y + 'n 
and r/o and made the identifications 

V'xO := V'2y+(0,r,) , V'xl := ^5y+(l,r,) 



(59) 



(60) 



in agreement with the relations (52) which show that when ?7o = 1 the operator 
ip2y+ri is a creation operator. 

In this notation the commutation relations (53) become 



{^l'r,"V^x,?o} = S^x'xf^r; 



Next we must determine a matrix N[ such that 



(61) 



64 $^(g-)y'(iV*)y'y (^+)y =EE^W (r!,,(iV,)y,yP+r,) ^2y+, (62) 

= J2Y1 ^2y'+r?' (^t )y'»?',y'7 ^2y+^ • (63) 



y'>y 
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In the above equation color taste and Dirac indices have been omitted. We observe 
that 



(7o7fc ® 1) Pt^^P+Tr, = ^Sor,, [(7o7fc ® 1) ± (7075 ® htk)] Tr, 



1 



1 T i-iT' 



r]+0-kJ 



and 
so that 



tr 
and 



r;,(7o7fc ® 1) Plt^P+'^v ='^SovoK'o^,k 



rj,(70 ® 1) P+r^ = 6or,oSlr,',^l,Tn 

iT(-l)^^ 



(64) 
(65) 

(66) 



2 i^OrikSri',ri+k + ^^Vk^r^' ,7]-%) 

(67) 



tr 



r;,(7o®i)p+r, 



4 I^Orjo'^liJo^r?'?? 



(6^ 



Finally we get the transformed iV-matrix 



f 1 ' 

Notice that the terms that involve the gauge variables refer to sites belonging to 
different blocks, while in the other terms the sites belong to the same blocks. The 
same operator can be re-labelled by using the coordinates x and r/o, then 



(AT'; 



\ 1 ' 



;(x')-5x',x+A<(x)] (70) 



where the values 77^ = 0, 1 simply control the presence of the gauge field according 
to the definition of u' given in (46). 
In conclusion 

q.Ntq+ = ^iN^i^o. (71) 

It should not be necessary to repeat that the expression of the transfer matrix so 
obtained is positive definite and performs time translations by two lattice spacings. 
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4.2 Coherent states 



In order to complete our analysis we perform the transformation also on the co- 
herent states. This will enable us to make, as a crosscheck, the derivation of the 
Lagrangian (47) starting from the transfer matrix. 
Let 



\a. 



(3) ■.= expl-2'J2J2K(^ 

I y 7:C 



(72) 



be a coherent state in the flavour basis, where a'^'^ and (3J'^ are Grassmann variables, 
such that 



Now 



2y+f? 



y 7:C y,»y 

y 7.C y,v 
and therefore, because of the anti-commutation relations (61) 

V^xo |a, /3) = E ^^y+'y ^Or,o |a, /?) = 2 tr (^rfp^^^tty j |a, /3) 
V'xi |a, /3) = E ^2y+r? '^i'?o |a, /3) = 2 tr (/3yr(i,^)) |a, /3) . 



This means that we can define 



ct^ :— 2 tr ( rlg^^-jay ■ , 



2 tr (/3yr(i,^)) 



and re-write 



|a, /3) = exp 



|0) 



(73) 

(74) 
(75) 

(76) 
(77) 

(78) 
(79) 



Notice that the Grassmann variables a, (3 and a' as well are defined at even times. 
The variable (3' instead, because of the matrix '^{i^r^) in its definition, must be 
considered attached at odd times. This is confirmed by the evaluation of the 
partition function using the transformed transfer matrix and coherent states. After 
the identifications 



^2x0+1 

we get the Lagrangian (47). 



(«2xo)* ' V^2xo = (/3Lo+l)' 
1^2x0+3 ' V'2xo+l = '^2a:o+2 



^0) 
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5 Conclusion 



Numerical simulations with Kogut-Susskind fermions are faster in the spin basis 
than in the flavor basis. Such calculations are usually performed in the lagrangian 
formulation, but we are interested in numerical simulations in the framework of 
the nilpotency expansion, that makes use of the transfer matrix. So we need an 
expression of the transfer matrix in the spin basis. In any case the knowledge of a 
positive definite transfer matrix in the spin basis is per se interesting being related 
to the unitarity of the theory. 

We found in the literature essentially two formulations of the transfer matrix 
in the spin basis. In the first one the Lagrangian is reduced by defining fermion 
fields and their conjugates at the odd, respectively even sites, and a transfer matrix 
is constructed that performs time translations by 2 lattice spacings [11, 12]. The 
fermion determinant even at vanishing chemical potential, is, however, not positive 
definite, which makes this way less suitable to numerical simulations. 

In the second formulation [11], a positive definite transfer matrix, called T^, 
was defined that also performs time translations by 2 lattice spacings. As a conse- 
quence the corresponding Fock space must be constructed on blocks. The explicit 
construction of such Fock space, however, is not given. 

If the Fock space is associated to a block, we can get the transfer matrix in 
the spin basis by a unitary transformation from that in the fiavor basis, whose 
expression, together with the construction of the Fock space, are known. The 
transfer matrix in the fiavor basis is expressed in terms of a matrix , and the 
transformed matrix is given in terms of the matrix A^', given explicitly in (70). 
In order to do numerical simulations in the nilpotency expansion all we need is 
to replace everywhere in the equations of the nilpotency expansion A^ by A^' and 
remember that the gauge fields are now defined on blocks. 

It would be now natural to compare our result with the expression of the 
previously derived transfer matrix [ ']. One might expect that such a comparison 
should provide the definition of the Fock space in the latter. Unfortunately this 
is not the case. The transfer matrix of [11] cannot be related to ours in a simple 
way, the most remarkable differences being that there is no requirement concerning 
the gauge variables which remain defined on the links of the original lattice, and 
creation and annihilation operators appear not only in exponential form but also 
as powers. 
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